New further integrability cases for the Riccati equation 
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New further integrability conditions of the Riccati equation dy/dx = a{x) + b{x)y + c{x)y^ are 
presented. The first case corresponds to fixed functional forms of the coefHcients a{x) and c{x) 
of the Riccati equation, and of the function F{x) = a{x) + [f{x) -— b^{x)] /Ac{x), where f{x) is an 
arbitrary function. The second integrability case is obtained for the "reduced" Riccati equation with 
b(x) = 0. If the coefficients a{x) and c{x) satisfy the condition ±d\/ f{x) /c{x) / dx = a{x) + f{x), 
where f{x) is an arbitrary function, then the general solution of the "reduced" Riccati equation can 
be obtained by quadratures. The applications of the integrability condition of the "reduced" Riccati 
equation for the integration of the Schrodinger and Navier-Stokes equations are briefly discussed. 
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INTRODUCTION 



Recently two integrability conditions of the Riccati equation 

^ = a{;x) + b{x)y + c{x)y^ , 
ax 



(1) 



where a, b, c are arbitrary real functions of x, with a,b,c G C°°{I), defined on a realinterval / C 3? 0, have been 
obtained in Q. By introducing a solution generating function f{x) G C°°{I), one can attach to Eq. (|T]) the auxiliary 
Riccati equation 



dx 



d 

dx 



2c{x) 



b\x) - fix) 
Ac{x) 



+ b{x)y + c{x)y^. 



(2) 



Then, if the coefficients a, b, c of the Riccati Eq. ([T]) satisfy the condition 



a{x) 



d_ 

dx 



^bjx) ± ./Jjx) 
2c{x) 



b\x)~f{x) 
4c(x) 



(3) 



the general solution of Eq. ^ is obtained as 



y±{ 



-b{x) ± y/fjx) 
2c{x) 



(4) 



where C± are arbitrary integration constants 0]. 

In 3 the general solution of the Riccati Eq. ([2]) was obtained and discussed in detail in two cases. First, the 
functional forms of the functions b{x), c{x) and f{x), was fixed, and the general solution of the Riccati equation was 
obtained with a{x) given by Eq. ([3]). In the second case, by fixing the functional forms of the functions a{x), b{x) and 
f{x), the general solution of the Riccati Eq. ^ was obtained with c{x) given by Eq. 

It is the purpose of the present paper to consider some further integrability conditions of the Riccati Eq. ([l} , by 
using the approach introduced in Q . An integrability condition can be obtained by fixing the functional forms of the 
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functions a{x), c{x), and of the function F{x) = a{x) + [f{x) — /Ac{x), which allows us to obtain the general 

solution of the Riccati Eq. with b{x) determined from Eq. ([3]). 

An integrability case for the "reduced" Riccati equation with b{x) = is also obtained. Some physical applications 
of the solution are considered for the case of the Schrodinger equation, and of the Navier-Stokes equation, respectively. 

The present paper is organized as follows. The general solution of the Riccati equation with given a, c and F is 
derived in SectionUll An integrability case for the "reduced" Riccati equation is presented in Section HIIl Some physical 
applications of the solution generating method for the "reduced" Riccati equation arc presented in Section ITVl We 
conclude our results in Section fVl 



II. GENERATING SOLUTIONS OF THE RICCATI EQUATION BY FIXING a{x), c{x) AND F{x) 
Eq. ([3]) can be rewritten as 



F{x) 



dx 



-b{x) ± 

2c{x) 



where we have introduced a new generating function F (x) defined as 

f{x)~b^{x) 



F{x) = a{x) + 



Ac{x) 



From Eq. (jH]) we obtain f{x) as 

f{x) = Ac{x) [F{x) - a{x)] + b'^{x). 
Integrating Eq. ([5]) yields the following equation: 

2c{x)(^j F{x')dx' + F^^b{x)^±^J{x), 

where Fq is an arbitrary integration constant. By inserting Eq. ([7]) into Eq. ([8]), we find 

F {x') dx' + Fa 



2c{x) 



b{x) = ± ^/4c{x) [F{x) - a{x)]+b^{x). 



By squaring both sides of Eq. ^ , we obtain the following equation: 



2c(x) / F {x') dx' + Fo + b{x) 



= 4:c{x) [F{x) - a{x)] + b^{x). 



Now rearranging the terms of Eq. ([TU]) , it is easy to show that 

F{x) - a{x) 



b{x) 



c{x) 



f F {x') dx' + Fa 
By substituting Eq. (fTTj) into Eq. ([8|), Eq. ([8|) can be expressed as 



± = <x) 



F (x) dx' + Fo 



F {x') dx' + Fa 
F{x) - a{x) 



J'' F{x')dx' + Fa 

By substituting Eq. (fTTj) into the Riccati Eq. ([l]), the Riccati Eq. ^ takes the form 

F{x) - a{x) 



dy , . 
— = a[x) 
dx 



f F {x') dx' + Fo 



c{x) 



F {x') dx' + Fo 



y + c{x)y 



(5) 

(6) 
(7) 
(8) 

(9) 

(10) 

(11) 
(12) 

(13) 



Therefore we obtain the following: 

Theorem. Let F{x) g C°°{I) be an arbitrary function defined on a real interval / C 5ft. If the coefficient b{x) of 
the Riccati Eq. ([1]) satisfies the condition (jlip . then the general solution of the Riccati Eq. (|13p is given by 



y{x) = 



Co- j c{x')e I. L 
where Cq is an arbitrary integration constant. 



1 



+ / F{x')dx' + Fo, 



(14) 



A. Example 1 



The coefficients of the Riccati equation 
where fci, fc2, m and n arc arbitrary constants, satisfy the condition (jlip . A particular solution of Eq. (jisp is 



yp = Fo 



1 + m 



The general solution of Eq. is given by 
= Fo + — + 



exp| 










(m+l)(m+n+2) ' n+1 


} 






exp|fc2?"+^ 


r fci5-+i 1 

(m+l)(m+Ti+2) 


Fo 
n+1 





B. Example 2 



The coefficients of the Riccati equation 

dy J 

— ~ kix + 
ax 



(n — ki) X 
nx^/2 + Fo 



- 1 



1 



71X2 /2 + Fo 



satisfy the condition pT|) . It has the particular solution 



The general solution of Eq. (fTS]) is given by 



2/(2;) 



(nx^ + 2Fo 



nx 



Co-2/"e?«2 + 2Fo)-"^/"d^ ' 2 
C. Example 3 



The coefficients of the Riccati equation 



dy_ 

dx 



kx" 



ps'-'+VCs + 1) + Fo X 



px" — kx^' 



X px-'+^/{s + l) + Fo 



px 



s+1 



S + 1 



Fo 



where fc, m, n, s and p are arbitrary constants, satisfy the condition (|lip . It has the particular solution 



yp^ Fo + 



px 



s+1 



S + 1 



The general solution of Eq. ([?!]) is given by 

y(x) = -G(x) 

where the function G {x) is defined as 



Co - G {x) ' 

Fo(s + 1) (fcx™+i + Fo + Fos) + " - 2F1 (1, ^^i±2±i; 



m+n+s+2 . _ pa: 



s+1 ' sFo+Fo 



p.j,s+l +Fo+ FqS 



m + n + 1 



where 2F1 (^l, ; '"+"^^^"+2 . jg the hypergeometric function 2F1 (a, 6; c; z) 
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III. THE "REDUCED" RICCATI EQUATION: THE CASE b{x) = 
If the function b{x) = 0, the Riccati Eq. ([!]) takes the "reduced" form 

- a(x) + c(x)y^ (25) 



dx 

We assume that a particular solution yp{x) of Eq. (|25l) satisfies the condition 

^=a{x)+f{x), (26) 

where f{x) G C^{I) is an arbitrary function. Substitution of Eq. (pS)) into Eq. fixes the particular solution yp{x) 
as 



..(.)^±^f|. (27) 

Therefore we obtain an integrability condition of the reduced Riccati equation, alternative to Eq. (jS]), and which 
can be formulated as 



Therefore we have obtained an integrability condition for the "reduced" Riccati equation expressed as the following 
Theorem. Let f{x) G C°°{I) be an arbitrary function defined on a real interval / C 3?. Then the general solution 
of the "reduced" Riccati equation 



^=±;;^J44-/(.) + c(x),^ (29) 
ax ax y cyx) 



is given by 



y±{x) = ±J^ + , , , (30) 



where C± are arbitrary integration constants. If c(x) = 1, the solution of Eq. (|53p can be written as 

d 



y±{x) = ±\/f{x) - — In 
dx 



C, - / e±2/^ Vf(^}'i-"dx' 



(31) 



IV. APPLICATIONS IN PHYSICS 

In the present Section we consider some physical applications of the obtained integrability conditions of the Riccati 
equation. In particular, we will consider the application of the obtained integrability conditions to the case of the 
Schrodinger - Riccati system, and of the Navier-Stokes equations. 

A. The Schrodinger-Riccati system 

The one-dimensional Schrodinger equation for a potential V{x), 

il;"{x) + [E-V{x)]i^{x) =0, (32) 

where E = constant is the energy, the Planck constant h ^ 1, and the mass m of the particle is normalized to m = 1/2, 
by means of the transformation 
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can be transformed to a "reduced" Riccati equation of the form Q 

u'{x) = E -V{x) +u^{x). (34) 
Therefore if the solution of the " reduced" Riccati equation (p4|) is known, the wave function is given by 



ip{x) = ipo exp 



u (x') dx' 



(35) 



where ipo is an arbitrary constant of integration. 

For the Schrodinger - Riccati equation p4p c{x) = 1 and a{x) ^ E — V{x). Therefore for any given arbitrary 
function f{x) the Schrodinger - Riccati equation has a closed form solution, with the potential fixed by 



V{x)^ET^Vm + f{^). 



(36) 



As a first application of the integrability case given by Eq. we consider the choice f{x) — /qx" for the function 
f{x), where /o and n are arbitrary constants. Then the potential for this problem is given by 



V{x)=E + fox"T^VTox"/''-\ 
and the general solution of the Schrodinger - Riccati Eq. ([M)) is 

n/2 , (n + 2) e-^ 



u±{x) = ±V/oa;" 



(37) 



(38) 



where En (z) = e ^*dt /t" is the exponential integral function [1| . The wave- function corresponding to this potential 
is given by 



V'i i^) = ^±0 exp [t2//^x"/2+V(" + 2)' 



C^ + 2xE^ U^x-/'+']/{n + 2) 



(39) 



where ip±o are arbitrary integration constants. As a particular example we explicitly present the solutions correspond- 
ing to the case n — 2, for which the quantum potential is 



Vix) = foX^ + ETVfo, 



(40) 



corresponding, from a physical point of view, to harmonic motion [6j. The solution of the Schrodinger - Riccati 
Eq. (p4|) for this potential is 



(41) 



where (z) = erf(iz)/i and F_ (z) = erf(z), with erf(z) = {2/y/n) e * [J]. The wave function for the harmonic 
double well potential is obtained as 



ip± (x) = ip±o exp 



T Vfox'/2 C± + V^F± ( /(j/'x ) /f^ 



pl/4 



(42) 



Eqs. (PT|) and (|^^ have been used in the physical context of quantum mechanics since 1984 0, 13- 

As a second example of the application of the integrability condition given by Eq. p8|) to the Schrodinger - Riccati 

system we consider the case f{x) = fox" — E, n ^ 0. The corresponding potential, satisfying the integrability 

condition of the Riccati equation, is 



V{x) = fox" T 



fon X 



n-l 



2 ^fox" - E 

The general solution of the Riccati equation Eq. (l34l) for this potential is 



u± 



(x) = ±Jfox'^ -E + , 

^ ^ C±- J H±{x')dx' 



(43) 



(44) 
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By substituting Eq. (1441) into Eq. (|35t . the wave function can be represented by 

i>±{x) = '4>±oH^ (x) 
where we have introduced the function H± (x), given by 



C±-J H±{x')dx' 



i?±(x) =exp|±2xyi^ 



E 



, 1 1 1 fox'- 
2 + n2Fi 1,- + -;! + -; ^ 
2 n n E 



/{2 + n) 



(45) 



(46) 



In Eq. (|43)) for ah n ^ the potentials are singular at the points fox^ = E. This singularity can be avoided, from 
a physical point of view, by taking into consideration only bound physical states, that is, those for which the energy 
E is negative, and by considering only half-line (radial) cases. 



B. The Navier-Stokes equation 

The Navier-Stokes equation for a steady viscous flow, 

p (w • V) tJ = — Vp + pf + /^Aw, 



(47) 



where p is the density of the fluid, p is the fluid pressure, v is the velocity, fi is the dynamic viscosity, and / is 
the external force acting on the fluid, is one of the most complex equations of mathematical physics. Therefore 
reducing it to a simpler form, or establishing, by using some physically reasonable assumptions, a connection between 
the Navier-Stokes equation and some other equations of the mathematical physics, is of fundamental importance in 
obtaining some exact solutions of the Navier-Stokes equation. By introducing a strain field 7 into the velocity field 
u of the fluid so that u = [— {j/2))x, — {j/2))y, (7) z + {—dtp/dy, dtp/dx^)], where tp, 7 and W are functions of x, 
y and t, but not of z, solutions of the Burgers type can be constructed [9[. If the strain rate 7 is a function of time 
only, 7 = j{t), then it is related to the pressure p of the fluid by a Riccati equation of the form Q 



J+7'+P..(t) = 0, 



(48) 



where Pzz{t) is the second partial z-derivative of the pressure, which must be spatially uniform, a constraint necessary 
for the existence of stretched vortex solutions of the 3D Navier Stokes equations with uni-directional vorticity, which 
is stretched by a strain field that is decoupled from them. Therefore, if there is a function f{t) so that the pressure 
can be represented as 



the strain rate ^{t) for a viscous fluid flow can be obtained as 



7±(t)=±v^+^ln 



(49) 



(50) 



where F-t are arbitrary constants of integration. 

In [U and ^ it was shown that along a streamhne the two-dimensional Navier-Stokes equation can be written 
as a Riccati equation of the form 



ill — a[s)u\ + /3{s) — 0, 



(51) 



where the dot denotes the derivative with respect to the parameter s, a — l/2iy, where i> — fi/p is the kinematic 
viscosity, and /3 = — (1/z/) {q/ p — fi) s + C/u, with C an arbitrary constant. The parametrization of the trajectory 
is of the form $ : s — > (c/)i(s), (/!)2(s)) = {x,y), and u — v o (j) and q = p o (f), respectively. A solution of the Riccati 
equation in terms of the Airy function was obtained in . By using the integrability condition Eq. ([28| more general 
solutions of the two-dimensional Navier-Stokes equation can also be constructed. Hence, if there exists a function 
/(s) so that the condition 



± 



ds y a{s) 



(52) 
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is satisfied for all a{s), /3{s) and f{s), the general solution of Eq. (|51|) is given by 



U±i{s)^±.U^ + , , , (53) 

\ a{s) u±- ra{s')e^^f' ^ f(s")^(s")ds" 



where U± arc arbitrary integration constants. 



CONCLUSIONS 



In the present paper, by extending the work initiated in Q, we have obtained two integrability conditions for the 
Riccati equation, one for the "full" equation, and one for its "reduced" form, respectively. Both integrability cases are 
based on the correspondence between the initial Riccati equation and a more general equation containing a solution 
generating function f{x). If the coefficients of the Riccati equations and the function f{x) satisfy some differential 
integrability conditions, the general solutions of the considered Riccati equations can be explicitly obtained. However, 
we would like to mention that generally the linear term in the Riccati equation, containing the function 6(a;), can be 
eliminated from the initial Riccati equation by means of the transformation 

y{x)=e^^^^''>'''v{x), (54) 
with v{x) satisfying the "reduced" Riccati equation 

dV , . - b(x')dx' I / N f'bix'^dx' 2 /rr\ 

— = a{x)e J \ ) + cixje-' ^ ' v . (55) 
ax 

The integrability conditions given by Eqs. (|29p and ((53)) can also be applied to Eq. (|55|) . However, sometimes in 
practical applications the analysis of the solutions of the Riccati equation can be done easier in the general form, 
especially in the situations in which the integral of the function b cannot be obtained in an exact analytical form. 

We have also presented some explicit physical applications of the integrability conditions for the Schrodinger- Riccati 
and the two-dimensional Navier-Stokes equations. Therefore, for the evolution equation of natural processes that can 
be reduced to a Riccati form, the presented integrability conditions open the possibility of finding explicit exact 
solutions for physical models showing a complex dynamical behavior. 



Acknowledgments 

We would like to thank to the anonymous referee, whose comments and suggestions helped us to significantly 
improve our manuscript. 



A. D. Polyanin, V. F. Zaitsev, Handbook of exact solutions for ordinary differential equations, Boca Raton, Chapman & 
Hall/CRC (2003). 

M. V. Soare, P. P. Teodorescu and I. Toma, Ordinary differential equations with applications to mechanics, Dordrecht, 
Springer (2007). 

M. K. Mak, T. Harko, New integrability case for the Riccati equation, Applied Mathematics and Computation 218, 10974- 
10981 (2012); arXiv: 1204.6546 (2012). 

M. Abramowitz, I. A. Stegun, Handbook of mathematical functions with formulas, graphs, and mathematical tables, Wash- 
ington, D.C., National Bureau of Standards (1972). 

S. Fraga, J. M. Garcia de la Vega, and E. S. Fraga, The Schrodinger and Riccati equations, Berlin; New York, Springer- 
Verlag (1999). 

E. Merzbacher, Quantum mechanics. New York, Wiley (1998), p. 149-159. 

B. Mielnik, Factorization method and new potentials with the oscillator spectrum, J. Math. Phys. 25, 3387-3389 (1984). 
H. C. Rosu and J. Socorro, One-parameter family of closed, radiation-filled Friedmann-Robertson-Walker quantum uni- 
verses, Phys. Lett. A 223, 28-30 (1996). 

J. D. Gibbon, A. S. Fokas, and C. R. Doering, Dynamically stretched vortices as solutions of the 3D Navier-Stokes equations, 
Physica D: Nonlinear Phenomena 132, 497-510 (1999). 

G. Argentini, Exact solution of a differential problem in analytical fluid dynamics: use of Airy's functions, 
|arXiv:math/0606723 (2006). 

V. Christianto and F. Smarandache, An Exact Mapping from Navier-Stokes Equation to Schrodinger Equation via Riccati 
Equation, Progress in Physics 1, 38-39 (2008). 



